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ABSTRACT 

We examine the local conditions for radiative damping and driving of short wave- 
length, propagating hydrodynamic and magnetohydrodynamic (MHD) waves in static, 
optically thick, stratified equilibria. We show that so-called strange modes in stellar 
oscillation theory and magnetic photon bubbles are intimately related and are both 
fundamentally driven by the background radiation flux acting on compressible waves. 
We identify the necessary criteria for unstable driving of these waves, and show that 
this driving can exist in both gas and radiation pressure dominated media, as well as 
pure Thomson scattering media in the MHD case. The equilibrium flux acting on opac- 
ity fluctuations can drive both hydrodynamic acoustic waves and magnetosonic waves 
unstable. In addition, magnetosonic waves can be driven unstable by a combination 
of the equilibrium flux acting on density fluctuations and changes in the background 
radiation pressure along fluid displacements. We briefly describe the conditions under 
which these instabilities might be manifested in both main sequence stellar envelopes 
and accretion disks. 

Subject headings: accretion, accretion disks — instabilities — MHD — stars: oscillations 

1. Introduction 

It is well known that the presence of a substantial equilibrium radiation pressure gradient 
can destabilize optically thick astrophysical media. Prendergast & Spiegel (1973) first speculated 
that compressible fluid flow subject to a large radiative flux would be unstable to the formation 
of buoyant rarefied regions of enhanced radiation pressure, better known as "photon bubbles." 
The existence of these instabilities has never been fully demonstrated in hydrodynamic systems, 
but radiation pressure acting in stellar envelopes with opacity peaks can drive certain oscillation 
modes (now known as "strange modes") unstable (e.g. Wood 1976; Saio, Wheeler, & Cox 1984; 
Gautschy 1993; Kiriakidis, Fricke, & Glatzel 1993). Strange modes have also been found to exist in 
hydrodynamic models of accretion disks (Glatzel & Mehren 1996, Mehren-Baehr & Glatzel 1999). 

Gautschy & Glatzel (1990) showed that strange mode instabilities satisfy a "non-adiabatic 
reversible" (NAR) approximation, in which the modes occur in a medium with effectively zero 
specific heat capacity, resulting in a vanishing luminosity perturbation. A physical model for 
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strange mode instabilities within the NAR approximation was established by Glatzel (1994). He 
noted that opacity peaks act as a reflection barrier for high overtone acoustic modes, sonically 
decoupling different radial regions in the star. The basic mechanism for the instability was the 
relative phase shift between the pressure and density perturbations, resulting from the constraint 
of having zero luminosity perturbation. Further aspects of the physics of strange modes have 
been elucidated by Papaloizou et al. (1997) and Saio, Baker, & Gautschy (1998). In numerical 
studies, Shaviv (2001) has found that both standing and propagating acoustic waves are unstable 
in atmospheres with sufficient radiation pressure support, even when the opacity is pure Thomson 
scattering. 

Radiation pressure driven instabilities occur in magnetized systems as well. Arons (1992) iden- 
tified and investigated photon bubble instabilities in strongly magnetized atmospheres of accreting 
X-ray pulsars. He showed that such instabilities are caused by the enhanced buoyancy that occurs 
when radiation diffuses into relatively rarefied regions in compressible perturbations. With accre- 
tion disk applications in mind, Gammie (1998) performed a linear, local stability analysis of a static, 
magnetized, stratified medium and found an instability over a finite range of wavenumbers that he 
identified as being similar in nature to photon bubbles. Using the simplifying assumption that gas 
and radiation were coupled together just by Thomson scattering, Blaes & Socrates (2001, hereafter 
BS01) found that Gammie's (1998) instability could be extended to arbitrarily short wavelengths 
where it manifested itself as an overstable slow mode. In addition, the fast magnetosonic mode 
could also be unstable for sufficiently strong radiation fluxes, albeit with smaller growth rate. 1 
Begelman (2001) has constructed a periodic shock train solution that may represent one possible 
nonlinear outcome of the slow mode instability. 

The purpose of this paper is to investigate the local, WKB version of these compressive in- 
stabilities in a unified manner. We focus on the local driving of propagating waves, and do not 
address here the issue of whether and how this local driving can manifest itself as a global insta- 
bility of a standing, normal mode of the object in question. We treat the thermodynamics of the 
medium quite generally, making no restrictive assumptions on the equality of the gas and radiation 
temperature in the perturbed state or using the NAR approximation. We find that hydrodynamic 
instabilities driven by radiative diffusion exist under a wide range of conditions, for both gas and 
radiation pressure dominated media. However, a medium supported purely by Thomson scattering 
opacity exhibits no local driving of acoustic wave instabilities, in contrast to the conclusions of 
Shaviv (2001). On the other hand, the anisotropic stress produced by magnetic tension widens the 
applicability of these instabilities to even broader classes of equilibria, even those which have opac- 
ities given by Thomson scattering only. The MHD instabilities can exist even when the equilibrium 
magnetic pressure is less than either the gas or radiation pressure. 



1 BS01 also found that Alfven waves could be destabilized, an effect that is almost certainly due to the fact that 
they considered a rotating equilibrium and these waves then have a small compressible component because rotation 
couples them to slow modes. 
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All these instabilities may grow on the dynamical time scale or even faster. The instability 
mechanism originates from the interaction of the equilibrium radiative flux with density and opacity 
fluctuations in the perturbed flow, as well as changes in radiation pressure along fluid displacements. 
The ultimate source of free energy is provided by the stratified radiation field. 

This paper is organized as follows. In section 2, we describe our basic equations and assump- 
tions, concentrating in particular on the thermodynamics of the coupled gas and radiation fluids, 
and also derive the general dispersion relation for waves in a stratified, magnetized medium within 
the local (WKB) limit. In section 3 we discuss the solutions of this dispersion relation in the hy- 
drodynamic limit where there is no magnetic field, generalizing and extending the work of previous 
authors. We discuss in some detail the physics underlying these instabilities in section 4. We then 
incorporate the effects of magnetic fields in section 5, and in section 6 discuss how the physics 
of magnetoacoustic wave instabilities is related to the hydrodynamic instabilities of the previous 
sections. In section 7 we briefly discuss astrophysical applications of these instabilities to accretion 
disks and stars, deferring more detailed applications to observed phenomena for future work. Sec- 
tion 8 summarizes our main conclusions. Our WKB analysis in the main body of the paper is not 
completely rigorous, and we provide more solid mathematical justification for our conclusions for 
two particular cases in an appendix. We also include an additional appendix where we examine the 
effects of radiative diffusion on the magnetorotational (MRI) instability in differentially rotating 
flows, generalizing our previous work (BSOI) to include the more generic thermodynamics we use 
here. Readers not interested in the details of the linear perturbation theory may wish to skip ahead 
to sections 7 and 8 directly. Readers interested in obtaining a basic physical understanding of the 
causes of these instabilities may wish to focus on sections 4 and 6. 



2. Equations and Assumptions 

The basic equations of radiation magnetohydrodynamics (RMHD) have been discussed by 
Stone, Mihalas, & Norman (1992), and we adopt these equations here with slight changes in notation 
as well as some changes in the physics. The fluid equations we consider are 

g + V(pv)=0, (1) 

p^+v-Vv) =-Vp + /9g+i-(VxB)xB + ^F, (2) 

du ,_, m4 / 4/cB^g hP \ 

— + v • Vu + 7«V • v = KjpcE - KppcaT* - K T pc ^ k E, (3) 

at V "i-cC m e c A J 

dE ^4^ „ . fAk B T R hP \ „ 

— + v ■ VE + -EV ■ v = -V • F — njpcE + K V pcaTt + k tP c — 5 E, 4 

at 3 V "ieC m e c z J 

0=-i VE -^F, (5) 
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and 

^- = Vx(vxB). (6) 

Here p, p, T g , and u are the density, pressure, temperature, and energy density in the gas, respec- 
tively. These quantities are related to each other by 



and 



U= —l (7) 

P=^, (8) 

where 7 is the ratio of specific heats in the gas and p is the mean molecular mass of the gas. 
We assume throughout this paper that 7 and p are constant, and thereby neglect the effects of 
composition gradients, ionization, and recombination. 

Other symbols in equations (l)-(8) have their usual meaning: c is the speed of light, a is the 
radiation density constant, k-Q is Boltzmann's constant, and m e is the electron mass. The vector 
v is the fluid velocity, B is the magnetic field, and g is the gravitational acceleration. We assume 
that g is time-independent throughout this paper, either because the gravitational field is due to 
some fixed external source (as is the case for a non-self-gravitating accretion disk), or because 
we choose to consider only short-wavelength perturbations that are well-described by the Cowling 
approximation . 

The radiation energy density E and the radiation flux F are defined in terms of frequency 
integrals over the radiation spectrum as measured in the local fluid rest frame: 

47r f°° 4ir 
E = — dvJ v = — J, (9) 

c Jo c 

where J u is the angle-averaged mean specific intensity in the local fluid rest frame, and 

> 

dvF v . (10) 



J 



We have simplified the RMHD equations by assuming that the radiation stress tensor is isotropic 
in the local fluid rest frame, i.e. that the tensor variable Eddington factor defined by Stone et 
al. (1992) is one third times the identity matrix. We have therefore neglected the effects of 
photon viscosity. We have also neglected terms corresponding to radiation inertia in the radiation 
momentum equation (5), leaving us with just the radiation diffusion equation. Such terms are 
generally negligible for the low frequency instabilities we explore in this paper (e.g. BS01). 

The two opacities kj and kf are also defined by integrals over the radiation spectrum in the 
local rest frame: 



and 



1 f°° 

^^— T du X ^( P ,T g )J u , (11) 
P J Jo 

1 f°° 

KF F = - / dv[ X f{p, T g ) + n c a T ]F u . (12) 
P Jo 
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Here x^iP^g) ^ the thermal absorption coefficient at frequency u, n e is the electron number 
density, and <tt is the Thomson cross-section. The Planck mean opacity is 

A _ /*0O 

KP = 'p~^J d "xf(p,T g )B u (T g ), (13) 
where B v {T g ) is the Planck function at the gas temperature. The Thomson opacity is 

- neCJT f1A\ 

k t = (14) 

and is constant due to our neglect of ionization/recombination processes and composition gradients. 
Note that we have extended the RMHD equations of Stone et al. (1992) to include electron scatter- 
ing. Coherent (Thomson) scattering merely allows the gas and radiation to exchange momentum, 
not energy. However, we have also allowed for incoherent (Compton) scattering in the last terms on 
the right hand side of the gas and radiation energy equations (3) and (4). We do this by treating 
Compton scattering in the frequency diffusion (Kompaneets) limit, neglecting anisotropies in the 
radiation field in the local fluid rest frame. A derivation of the Compton scattering terms can be 
found, e.g., in section 3.1 of Hubeny et al. (2001). The quantity u is an average frequency of the 
radiation field in the local rest frame of the fluid, defined by 



I J r 2 



The second term in parentheses represents the effects of stimulated scattering, and guarantees 
that there is no net heat exchange between the gas and the radiation when thermal equilibrium is 
established with the radiation field being blackbody at the local gas temperature. 

In principle, evaluation of the frequency integrals in the opacities and v requires a full solution 
of the radiative transfer equation coupled to the fluid flow. Rather than do that here, we adopt 
some simplified prescriptions that we believe will not alter the essential physics we wish to explore. 
First, we adopt the Eddington closure relation between the zeroth and second moments of the 
comoving frame specific intensity, consistent with our earlier assumption that the tensor variable 
Eddington factor is one third times the identity matrix. Second, because we will only consider 
perturbations about equilibria that are in local thermodynamic equilibrium (LTE), we assume 
that non-LTE departures in the perturbed flow can be adequately parameterized by a blackbody 
radiation spectrum at a temperature T r , which may be different from the gas temperature T g ? 

Under these assumptions, equation (15) implies that 

» = —. (16) 



2 Departures from LTE will of course generally not maintain a blackbody radiation field in the perturbed flow. 
Strictly speaking, we are defining T T as the effective temperature of the radiation field, and we are assuming that the 
various frequency moments we need are related to each other in the perturbed flow in approximately the same way 
as they would be if the radiation field were blackbody. 
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Equation (11) gives an expression for kj which more closely resembles the Planck mean opacity Kp, 



4vr 



dvxf(p,T g )B u (T r ] 



(17) 



pacTf J 

In addition, equation (12) can be replaced with an expression analogous to the Rosseland mean 

-l 



opacity, 



dv 



xl h (p,T g ) 



+ KT 



dB u 
dT r 



(18) 



Unlike kp, which is a function only of the gas density and temperature, kj and kf are also functions 
of the radiation temperature. Finally, the blackbody assumption implies that T r can be related to 
the radiation energy density through the equation of state 



E = aTt 



(19) 



To summarize, equations (l)-(8), (13), and (16)-(19) are the basic dynamical equations we will 
use throughout the rest of this paper. 



2.1. Equilibrium 

For the majority of this paper we will consider the dynamical stability of static equilibria in 
which the equilibrium fluid velocity is zero. As mentioned in the previous section, we also assume 
that the equilibrium is in LTE, i.e. that the unperturbed gas and radiation temperatures are 
identical, 

T g = T r = T. (20) 

This can be somewhat problematic for applications to accretion disk models, which in some cases 
can be effectively thin. However, we do not presume to model the turbulent dissipation that must 
exist in these flows. In fact, our purpose in that application is to consider equilibria that might be 
useful starting points for exploring the development of this turbulence by simulation. Lacking an 
explicit treatment of dissipation, we are forced to consider equilibria that are in LTE and in local 
radiative equilibrium. 

Under these assumptions, the only nontrivial partial differential equations describing our equi- 
librium are those expressing hydrostatic equilibrium, 

= pg- Vp+— (V xB)xB + — F, (21) 

47T C 

radiative equilibrium, 

= V-F, (22) 

and radiative diffusion, 

= -I V £-^F. (23) 

o C 
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Throughout this paper we will assume that the equilibrium magnetic field is uniform, so that the 
Lorentz force term in equation (21) vanishes. Equation (21) then has the immediate consequence 
that all thermodynamic variables are constant on horizontal surfaces, i.e. those surfaces that are 
perpendicular to the gravitational acceleration g. We define z to be a vertical coordinate that 
increases upward, and z to be the corresponding upward unit vector. Then g = — gz, with g > 0. 



2.2. The Linear Perturbation in Total Pressure 

When we consider the behavior of linear perturbations about the equilibrium just discussed, 
the problem naturally divides into two parts: the thermodynamics of the perturbed flow, and the 
coupling of this with the perturbed magnetic and velocity fields. We examine the former problem 
first by determining the perturbed total (gas plus radiation) pressure. 

We begin by considering expressions for the perturbations in kj and Kp. The functional 
dependence of these two opacities on p, T g and T r immediately implies that 

Skj dlnKjSp dlriKj5T g dlnKj8T T 

~k7 ~ dlnp ~p~ + dlnTg T g + 91nT r T r ^ ' 

and 

Skp d In up bp d In Kp 5T g 

~kV = dlnp ~p~ + dlnTg T g ' ^ ' 

(We use Eulerian perturbations throughout this paper, i.e. 5Q is the Eulerian perturbation in the 
quantity Q.) Here all opacity derivatives are evaluated in the equilibrium state which is in LTE 
with T T = T g and Kp = kj. Because of this, equations (13) and (17) immediately imply 

d In Kp d In kj 
dlnp dlnp 

and 



(26) 



din kp = d In k 3 din k 3 
<91nTg ~ <91nTg <91nT r ' 1 ' 

These two conditions on the absorption opacity derivatives greatly simplify the mathematical 
form of the thermal coupling between the gas and radiation in the perturbed flow. Linearizing the 
gas energy equation (3) and the radiation energy equation (4), and using equations (16), (19), and 
(26)-(27), we obtain 

dSu r e 4(5T T -5T R ) , x 

— + <5v • Vu + 7 uV • 5v = co a E y T ^ (28) 

and 

dt 3 T v ' 

Here cj a is an angular frequency describing the rate at which radiation and matter are thermally 

coupled in the perturbed flow, defined by 



w a = 



. , l91n Kj \ k B T 

Kp 1 + - m + K T- 



AdlnTrJ m e c 2 



pc = K a pc. (30) 
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Equations (28) and (29) can be further manipulated in a physically revealing manner. We 
eliminate V • <5v from them using the perturbed continuity equation (1), 

^ + pV • <5v + 5v ■ Vp = 0. (31) 

In addition, we eliminate SF using the perturbed radiative diffusion equation (5), which can be 
written as 

bp Jkf" 



5F = V5E - F — + — . (32) 

3kfP V P k f J 

Eliminating Su using equation (7), equation (28) becomes 

f - cf-^ = 4( 7 - - (7 - Dp«v ■ VS, (33) 

where c g is the adiabatic sound speed in the gas, 

1/2 



TpY 

p J 



(34) 



and S g is the entropy per unit mass in the gas, 

fee 

Scr = — r ln(pp~ 7 ) + constant. (35) 

g P{l ~ 1) 

Equation (29) can be written in a very similar fashion, 

9 A.rA 2 9Sp 4 (5T r -5T„) 1 mr / c 1 rrn \ 1^ /<5p 5k f \ 

— ( -^J-c^^f = --uj a E y r ^ — iZ__ p T(5v-V5 r +V-( ^— V-<5^J+-F-V f + — ^ J , 



^ V3 / r dt 3 T 3 r V 3k fP 

where c r is the adiabatic sound speed in the radiation, 



(36) 



and 5 r is the entropy per unit mass in the radiation, 

AE 

S T = — — + constant. (38) 
3pT v ; 

Note that we have simplified the terms involving the equilibrium flux F in equation (36) by using 
the radiative equilibrium equation (22). 

Thus far, we have not made any assumptions about the wavelengths of the perturbations we are 
considering, and we could proceed from this point by doing a full global perturbation analysis. We 
wish to focus on local instabilities in this paper, however. Hence from now on, we invoke the WKB 
ansatz by adopting a plane wave spacetime dependence oc exp[i(k • r — cut)] for all perturbations. 
Here r is the position vector of the point of interest, k is the perturbation wave vector, and uj is 
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the perturbation angular frequency. As a result, the time derivative operator d/dt in equations 
(33) and (36) is replaced by — iuj, and the spatial derivative operator V in the last two terms of 
equation (36) is replaced by ik, i.e. 



6Kpp 3 



+ 3 FV 



bp 8np 
P hf 



ck 2 
3kfP 



5E_ 
IT 



+ -ik • F 

3 



bp bnp 
P n F 



(39) 



The diffusive instabilities in which we are interested are driven by background gradients, in par- 
ticular the presence of a nonzero equilibrium radiative flux F = — (c/3kfp)VE. Hence this last 
step may at first sight appear to be dangerous: gradients in c/{3kyp) arising from the first term in 
expression (39) might be as important as the second term that depends on F as we take the short 
wavelength limit k — > oo. As may be verified a posteriori, however, it turns out that the insta- 
bilities have SE scaling as k _1 bp at high k, so that our ordering is in fact consistent. This makes 
physical sense, because at short wavelengths radiative diffusion is very fast and perturbations in 
the radiation temperature are smoothed out. 



After some algebra, equations (8), (19), (33), and (36) can be used to derive an equation for the 
total (gas plus radiation) pressure perturbation in terms of the density and velocity perturbations, 



Here 



bp + \bE = (Ac 2 g + Be 2 + C)bp - 



( 7 - l)AVS z + -BVS r 



A = 



UJ 

V 



UJ + 



ick 2 
3kfp 



+ iu a 1 + 



AE\ 
3p~) 



V 



k-F 



/e Tr 




V AE 


3p J\ 



4E\ 

uj + icj a (7 - 1) ( 3 + — J 



(40) 



(41) 
(42) 



C 



1 \ U^E 
V\ p 



4\ ick 2 . _ 
7 -3j W+ 3^ (7 - 1) 



k-F 



- (i + e p - e Tg ) + i( 7 - iK (i + @p - e rg - e Tr 



+i ( 7 - IK ( ^) (i + %) 



(43) 



and 



ick 2 k-F^ 

V = UJ + + —^QTr 

3kfp AE 



uj + 



4i( 7 - 1)uj & E ' 
p 



+ ioJ a 



k • F 

UJ + (7 - 1)— — 0Tg 



(44) 



The quantities are defined as logarithmic derivatives of the flux mean opacity with respect to 
the variable in the subscript, i.e. 



_ 8 In k f 



-, and &Tr = 



dlnT r 



(45) 



Equation (40) expresses all the thermal physics and is the one we shall use throughout the rest 
of the paper. 
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2.3. Coupling to Density, Velocity, and Magnetic Field Perturbations 

Now that we have determined the total pressure perturbation, the only work remaining is to 
couple this to the perturbed continuity, gas momentum, and flux- freezing equations. Employing 
the WKB ansatz, the perturbed continuity equation (31) becomes 

-iu5p + ipk ■ <5v + <5v • Vp = 0. (46) 

The perturbed gas momentum equation and flux freezing equations may be written as 

-iojp5v = -ik [5p + \dE ) + g5p + -^-(k x SB) x B (47) 

\ 3 / 47T 

and 

-iuSB = ikx (<5v x B), (48) 

respectively. 



2.4. The Dispersion Relation 

After some algebra, equations (40) and (46)-(48) may be combined to give a dispersion relation 
for short wavelength modes on a static, stratified and magnetized equilibrium: 



+ Co 2 (k 2 -k 2 z )(AclN 2 g +Bc 2 r N?) 
+ iuj 2 [oj 2 k - k 2 (k ■ v A )v A ] 
- [uj^ -u 2 k 2 v\-2uj 2 k z v Az k-v A + L0 2 k 2 z v% + k 2 {k-v A ) 2 v\ z ] g- Vlnp, (49) 

where 

C? = to 2 - (k- v A ) 2 , (50) 
v A = B/(47rp) 1 / 2 is the vector Alfven speed, 



-Vp(A-l) + ^-VE{B- 1) +CVln/> 
P 3p 



N s=S- (^ 1 Vp-Vlnp)=- il l)pT g-VS g (51) 
is the Brunt- Vaisala frequency in the gas, 

N?s - g -(^ VE - v ^=-w s ' vs ' (52> 

is the Brunt- Vaisala frequency in the radiation, and the quantities A, B, and C are defined by 
equations (41)-(44). 
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Equation (49) is a dispersion relation for eight modes. We originally had nine first order, 
time-dependent perturbuation equations (one continuity, three momentum, two energy, and three 
flux- freezing) , and so we expect nine modes. However, one of these modes has zero frequency and 
is inconsistent with the additional constraint of Gauss' Law that V • <5B = 0. In sections 3 and 5, 
we provide a full discussion of the unstable waves contained in equation (49). 



3. Hydrodynamic Instabilities 



Before considering the full effects of MHD, it is useful to first explore radiative diffusion insta- 
bilities in the hydrodynamic case. Setting va = in the general dispersion relation (49), we obtain 
the following equation for hydrodynamic modes: 







= to 4 - co 2 k 2 (Ac 2 R + Be 2 + C) + (k 2 
+ iuj 2 k ■ 



k 2 ) {Ac 2 N 2 + Bc 2 N 2 ) 



-Vp{A-l) + ^VE {B 
VP 3p 



1) +CVlnp 



us 2 g ■ Vlnp. 



(53) 



This equation can be written as a sixth order polynomial in u>, which is to be expected given that 
it arises from six time-dependent hydrodynamic perturbation equations. However, one power of oo 
can be factored out, indicating that one of the modes always has zero frequency. This mode has 
zero total pressure perturbation Sp + SE/3, at least when the wave vector is not entirely vertical. 



3.1. Short Wavelength Limit 



In the short wavelength limit, the remaining five modes described by equation (53) can be 
easily factored: 



2 4w a £(7-l) 
LU H OJ — 

IP 



1 



kz 



Nt 



(54) 



From left to right, the three factors respectively correspond to adiabatic acoustic waves in the gas, 
a purely damped radiative diffusion mode, and gravity waves in the gas modified by damping of 
gas temperature fluctuations by radiative emission and absorption. Note that radiation pressure 
and radiative buoyancy have been lost due to the rapid radiative diffusion at short wavelengths. 

The behavior and character of the gravity waves depends on the relative magnitude of the 
Brunt- Vaisala frequency and a characteristic gravity wave thermal coupling frequency 



LO thg 



2o; a £( 7 - 1) 
IP 



Solving equation (54) for the gravity mode frequencies, 



-iu thg ± {-uj 2 hg + N 2 



(k • z) 2 ] } 



1/2 



(55) 



(56) 
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If iV|[l - (k • z) 2 ] > w 2 hg , then the gravity waves are damped. If < iV 2 [l - (k • z) 2 ] < cj 2 hg , then 
this damping is so strong that the modes lose their wavelike character and have purely negative 
imaginary frequencies. If iV 2 < 0, then one of the gravity modes becomes convectively unstable, 
but with a growth rate that is reduced compared to ideal hydrodynamic convection. This reduction 
is due to the fact that emission and absorption damp gas temperature fluctuations in the wave, 
relative to the nearly uniform radiation temperature. 

The acoustic waves are more interesting. Expanding equation (53) to first order about the 
high k limit, we find 



The first two terms within curly brackets represent damping by radiative diffusion and emis- 
sion/absorption, respectively. The third term, which depends on the equilibrium radiative flux 
F, can give rise to instability if it dominates the first two. For a Kramers type opacity law in a gas 
with 7 = 5/3, p + (7 — l)0Tg < 0, and equation (57) then implies that downward propagating 
sound waves are potentially unstable, while upward propagating waves are damped. On the other 
hand, if there is no opacity perturbation, Q p = ®Tg = 0, and there is no acoustic instability. Such 
is the case for a medium with pure Thomson scattering opacity. 

Instability requires that the driving term dominate the two damping terms in equation (57), 
so that a rough, order of magnitude instability criterion is 



where = |0 P + (7 — l)©r g |. The first part of this criterion has an intuitive interpretation: if is 
of order unity, short wavelength acoustic waves will be unstable if the radiative flux is transporting 
the local radiation energy density faster than the sound speed in the gas. 

As we will see in section 3.3 below, the maximum growth rate of the instability, when it 
exists, occurs for wavenumbers such that the first term in equation (57) exceeds the second, i.e. 
for k > kf-F0/(cc 2 ). In a radiation pressure dominated medium, this implies that k > Og/c 2 , i.e. 
wavelengths shorter than the gas scale height will have maximal growth rates. This growth rate 
itself will be ~ @(g/c g ), faster than the reciprocal of the local free fall time by the ratio 0c r /c g . 

On the other hand, in a gas pressure dominated medium, the growth rate is ~ QE/(HTpc g ), 
where Ht is the temperature scale height. Taking this to be comparable to the pressure scale height 
~ c 2 /<7, we find a much smaller growth rate ~ Q(g/c g )(E/p), slower than the reciprocal of the local 
free fall time by the ratio 0(c r /c g ) 2 . This growth rate occurs for wavenumbers > Q(g/c^)(E/p), a 
threshold which is in violation of our WKB requirement that k S> g/c g for gas pressure dominated 
media. If the gas and radiation temperatures are not tightly locked, it appears difficult to produce 
vigorous acoustic wave instabilities in gas pressure dominated media by this mechanism. 



cj = ±kc e — 



. kf r f£cg 4£ccj a ( 7 - l) 2 



(k.F)[0 p + (7-l)0Tg]| + O(^ 1 ). (57) 



2cc g \ 3 Kppcg 




(58) 
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An expression similar to equation (57) appears to have been first derived by Hearn (1972; his 
equations 25, 26 and 30 and surrounding discussion), the only difference being that he does not 
have the radiative diffusion damping term [the first term in curly brackets in equation (57)]. This 
is because Hearn (1972) was interested in radiative amplification of acoustic waves in optically thin 
media, and therefore did not include the dynamical evolution of the radiation field in his analysis. 3 



3.2. Short Wavelength Limit With T g = T r 

If we first take the u; a — > oo limit in equation (53), so that the gas and radiation are tightly 
thermally coupled, then a mode is eliminated from the dispersion relation which is now only fifth 
order, including the zero frequency mode. The fact that we have lost a mode makes sense, as this 
limit corresponds to replacing the two time dependent gas and radiation energy equations with a 
total energy equation and the time-independent condition that the gas and radiation temperatures 
be equal. On taking the short wavelength limit, the four modes with nonzero frequencies factor as 
follows: 
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where 



Q = - 



1/2 



7 



1/2 



(59) 



(60) 



is the isothermal sound speed in the gas. 



These short wavelength modes resemble those of equation (54). Acoustic waves propagate at 
the isothermal sound speed, because the gas temperature is locked to a radiation temperature which 
is made nearly uniform by the rapid radiative diffusion. The two gravity modes have collapsed to 
a single mode in the last factor of equation (59). This mode has a purely imaginary frequency, and 
has therefore lost its wavelike character. It is unstable if 



^ N* + 
4(7 -1)£ g ^ 



1 9 

-iV r 2 < 0. 



(61) 



However, radiative diffusion strongly diminishes the growth rate of this instability in the short 
wavelength limit, with \u\ oc k~ 2 as k — > 00. This is because buoyancy is strongly suppressed: a 
perturbed parcel of gas always has the same density as its surroundings, because it has the same 



3 It is worth noting that Hearn's (1972) analysis had an error in his treatment of the gas energy equation (his 
equation 8). In this equation the partial time derivatives should be Lagrangian derivatives. It turns out, however, 
that his results are correct because he also neglected gas pressure gradient contributions to the equilibrium hydrostatic 
balance. These two errors cancelled one another. 
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pressure and is forced to have the same temperature due to the rapid radiative diffusion and rapid 
thermal coupling with the radiation. 

The acoustic waves are again unstable when first order corrections to their frequencies are 
made: 

" 'AE 



co = ±kc[ — i 1 + — — 

2cci V 4£ 



+ pj c iT (k-F) e, 



+ 0(k- 1 ). (62) 



The first term inside square brackets again represents damping by radiative diffusion, while the 
second term will drive instability if it is larger. For, e.g., Kramers law type opacities, 6 P > and 
now the upward propagating wave is unstable while the downward propagating wave is damped. 
Once again, a pure Thomson scattering medium exhibits no acoustic wave instability. 

An order of magnitude instability criterion from equation (62) is 

F@ p > m&x[E,p]ci, (63) 

which should be contrasted with equation (58). Note that there is no threshold from emission and 
absorption because the gas and radiation temperatures are equal. Here, acoustic waves are unstable 
if the local radiative flux is transporting the local thermal energy density, whether it be dominated 
by gas or radiation, faster than the gas sound speed. 

In a radiation pressure dominated equilibrium, equation (62) implies that the growth rate of 
the acoustic wave instability is ~ kfFQ p / \cci) ~ Q p g/c\ from hydrostatic equilibrium. In a gas 
pressure dominated medium, the growth rate is ~ n^F pC{Q p /{cE) ~ c\Q p /Ht- Hence the growth 
rate in a gas pressure dominated medium is also ~ Q p g jc\. In contrast to the previous case where 
gas and radiation did not exchange heat rapidly, radiation pressure support is not required to 
obtain high instability growth rates when the gas and radiation temperatures are the same. In 
a gas pressure dominated equilibrium, the small radiation pressure fluctuations produced by the 
damping and driving forces produce gas pressure fluctuations that are ~ (p/E) times larger just 
by the fact that the gas and radiation temperatures are locked together by rapid absorption and 
emission. 

Hearn (1972; equations 22 and 23) was also the first to derive an expression similar to equation 
(62), although again without the first damping term because he was interested in optically thin 
media. In addition, Hearn's (1972) approximations resulted in a replacement of the factor (1 + 
0.75p/E) multiplying the damping and growth rates in equation (62) with unity. This factor 
becomes very important in gas pressure dominated media. 

We note that our result that pure Thomson scattering media exhibit no local hydrodynamic 
acoustic wave instability disagrees with the conclusion of Shaviv (2001), who claimed that his "Type 
II" instability represents just such a local instability. In fact, this instability appears to have rather 
long vertical wavelengths. The growth rates actually depend on the boundary conditions , and also 
the location of the boundary (Glatzel 2003, private communication), indicating that it is global in 
nature, not local. 
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3.3. The Limit of Negligible Gas Pressure 

A number of authors have studied versions of these instabilities in the limit where gas pressure 
is completely negligible compared to radiation pressure. The gas sound speed then vanishes, and 
one then loses the acoustic wave nature of the instability in the short wavelength limit. 

Setting p = = c 2 in the hydrodynamic dispersion relation (53), and then adopting the ansatz 
that u oc k 1 ! 2 as k — > oo (Glatzel 1994, Gammie 1998), it is straightforward to show that there are 
two modes with 

w 2 = . FQp = _ ik . gQp ( 64 ) 

Assuming Q p > 0, one of these modes corresponds to an upward propagating unstable wave, while 
the other corresponds to a downward propagating wave which is damped. It is interesting to note 
that this result agrees with the thermally locked acoustic wave frequency, equation (62), if one 
squares that frequency and then takes the limit of negligible gas pressure. Note that the damping 
terms, which set the threshold of instability, vanish in the limit of negligible gas pressure. 

If we first assume negligible thermal coupling between the gas and the radiation (w a — > 0) , and 
then consider the zero gas pressure limit, we get instead two modes that reflect the form of the 
driving in the two temperature regime: 

J = • F [Q p + ( 7 - l)0 T g] = -ik ■ g [@ p + ( 7 - l)6 Tg ] . (65) 

Once again, the growth rate exhibits a k 1 ^ 2 dependence. Note that the damping terms again vanish 
in this particular limit. 

The acoustic wave instabilities we have been discussing throughout this section are the local, 
WKB versions of the strange mode instability discussed in the stellar oscillation literature. Glatzel 
(1994) has discussed a physical origin of strange modes in which he presents a WKB analysis of 
the growth rates of purely radial modes (i.e. vertical modes in our geometry) in the limit of zero 
gas pressure. His equation (5.8) is (nearly 4 ) identical to our equation (64). He did not recover 
the damping terms of equation (62) for two reasons: he assumed negligible gas pressure, and he 
also invoked the NAR approximation, a point that we discuss further below in section 4.1. Glatzel 
(1994) also presented a global analytic solution for unstable strange modes in which he showed 
that unstable and damped waves propagate in opposite directions (see discussion after his equation 
5.17), in agreement with our conclusions. 



4 Glatzel's equation 5.8 is lj 2 — ±i2gQ p k, and we actually only recover the version of this equation with a plus 
sign. In addition, we do not have his factor of two. The reason for both of these facts is that his definition of the 
wavenumber k differs from ours. Glatzel performed a WKB analysis on a variable IT which was a transformation of 
the Eulerian pressure perturbation (his equation 5.4). Our analysis only recovers modes with short wavelengths in 
the pressure perturbation, and this is true of only one sign of his dispersion relation. Moreover, our wavenumber for 
these modes is twice his wavenumber. 
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3.4. Numerical Results 

Table 1 summarizes the basic conclusions of the analytic work on acoustic wave instabilities 
presented in the previous two subsections. We now discuss numerical solutions to the full dispersion 
relation (53) which show how these instabilities operate in different regimes, illustrating the basic 
formulas in Table 1. 

At sufficiently short wavelengths, photons diffuse across a wavelength faster than the wave 
period, and radiation temperature fluctuations therefore become small. Whether or not the gas 
and radiation temperatures are locked together is therefore determined by whether absorption and 
emission are rapid enough to drive the gas temperature fluctuations to be small as well. From 
the perturbed gas energy equation (33), it is clear that the gas temperature will be locked to the 
radiation temperature provided the wave frequency lo is much less than a characteristic thermal 
frequency 

"4( 7 - 1)E~ 



P 



u a . (66) 



Figure 1 depicts unstable wave growth rates as a function of wavenumber for a radiation 
pressure dominated equilibrium, for different, small values of Uth- Although such equilibria are 
generally dominated by Thomson scattering opacity, for purposes of illustrating the physics, we 
have assumed Kramers type values for the flux mean opacity derivatives: © p = 1 and @Tg = —3.5. 
In any real application to a radiation pressure dominated equilibrium, these derivatives will be 
reduced approximately by the ratio of the absorption opacity K a = uj a /(pc) to kt, and this would 
lead to a reduction in the unstable growth rates, or possibly stabilization. We present results 
appropriate for more realistic applications later on below. 

At the low values of Uth shown in Figure 1, gas and radiation temperatures are decoupled at 
high wavenumber, and the asymptotic growth rates agree with equation (57). These unstable waves 
propagate downward because P + (7 — l)0Tg < 0. As wth increases, their asymptotic growth rates 
decrease and eventually damp. At the same time, a new set of modes appear at low wavenumber 
with growth rates that climb with increasing thermal coupling between the gas and radiation, as 
shown more clearly in Figure 1 (b) . These modes propagate upward and at higher thermal coupling 
turn into the one-temperature acoustic waves described by equation (62) at high wavenumber. 

This is shown in Figure 2(a), which depicts the growth rates as a function of wavenumber 
for high values of u; t h, for the same equilibrium parameters used in Figure 1. In contrast to the 
two-temperature behavior, the asymptotic growth rates at high wavenumber must eventually fail 
because at short enough wavelengths the modes must return to being adiabatic in the gas alone, 
and the assumption of thermal locking of the gas and radiation underlying equation (62) must 
break down. Figure 2(b) illustrates that for wavenumbers k > u; t h/c g , the real part of the phase 
velocity returns to the adiabatic gas sound speed. 



The one temperature instability growth rate cuts off at even lower wavenumbers than this, 
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however, and the reason appears to be a breakdown of thermal locking related to an interplay 
between the equilibrium radiation flux, absorption and emission, and fluid inertia. Numerical 
exploration shows that the growth rate cuts off for wavenumbers k J> (wth<7|l + ©p|/ c f) 1,/2 f° r 
both gas and radiation pressure dominated equilibria, and this formula for the maximum cutoff 
wavenumber is included in Table 1. 

For demonstration purposes, we have been considering equilibria with rather artificial opaci- 
ties. Figure 3 shows the unstable wave growth rate as a function of wavenumber for a gas-pressure 
dominated equilibrium with a realistic Kramers type opacity law, appropriate for mid- to up- 
per main sequence stars. For sufficiently high thermal coupling, as is indeed present in these 
stars (see section 7.2 below), unstable upward propagating waves exist over a broad range of high 
wavenumbers. Interestingly, long wavelength instabilities also appear in downward propagating 
waves. These downward wave instabilities exist because of the temperature dependence of the 
opacity (&Tg = —3.5). Because the gas and radiation temperatures are tightly locked, temperature 
fluctuations in the gas are smoothed out by radiative diffusion. Hence opacity fluctuations caused 
by 8T are smaller than opacity fluctuations caused by 5p. Therefore the temperature dependence 
of the opacity does not affect the growth rate at short wavelengths. This agrees with equation 
(62). However, at long wavelengths where radiative diffusion is slow, significant gas temperature 
fluctuations exist in the wave, and the resulting opacity fluctuations reverse the radiative driving. 
Note that these long wavelength, downward unstable waves do not appear to exist in radiation 
pressure dominated equilibria, as shown in Figure 2. 



4. Physics of the Acoustic Wave Instabilities 

The results of the previous section are based on a brute force analytic and numerical attack on 
the full dispersion relation. In this section, we present a more physically-motivated derivation of 
some of these results. (This turns out to be a useful starting point for understanding the unstable 
MHD waves that we discuss below in sections 5 and 6.) It is the coupling of the equilibrium radiation 
flux to opacity fluctuations that drives acoustic waves unstable. To see how this works physically, 
consider the forces acting on a perturbed parcel of gas in the acoustic wave. [This approach is 
somewhat analogous to using a zeroth order eigenfunction to compute damping or driving rates 
through a work integral in stellar oscillation theory (Cox 1980).] 

The gas and radiation energy equations (33) and (36) may be written in the WKB limit as 
tp = ff, - lfv ■ (Vp - 4V P ) + 4 ' h -^" E (gLlgk) (67) 

and 

r y_E _ c> _^ Sv . VSi _^e(st,-sta 

\ 3kfpuj J 3 3a; 6uj \ T J 
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(i + e P )| + e Tr f + e T! ^ 



(68) 



respectively. The radiation energy equation (68) immediately implies that for short wavelength 
(k — > oo) acoustic waves, SE ~ 0(k~ l )bp. Physically, at high k the radiation diffusion time 
across a wavelength is shorter than a wave period, so that radiation temperature fluctuations are 
smoothed out. On the other hand, because acoustic waves in the gas are adiabatic in this limit, the 
gas temperature fluctuation is given by bT g /T = (7 — l)6p/p. Hence the gas pressure perturbation 
is given by 

bp = clSp - % -b, ■ (V P - clVp) - 4 * (7 " l? ^ E ( ^ 

6 to 6 ' UJ \ p 

The first term on the right hand side is the dominant term at short wavelengths, and represents the 
acoustic response of the gas to compressive perturbations. The second term arises from buoyancy. 
The last term is a damping term caused by radiative emission/absorption: temperature fluctuations 
in the gas emit and absorb photons to try and come into equilibrium with the nearly uniform 
radiation temperature. 



+ 0(k 
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(69) 



To the same order in A; 1 , the radiation pressure perturbation is 
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bp + H (k • F) [1 + e p + (7 - l)6 Tg ] bp + o(k- 



(70) 



The first term on the right hand side represents damping by radiative diffusion: a density fluctuation 
will try to compress photons, but these diffuse quickly and transport heat out of the fluctuation. 
The last term is the most interesting. It is an extra radiation pressure force arising from the 
equilibrium radiation flux acting on density and opacity fluctuations in the gas. As we shall see, it 
is this extra force that is destabilizing if it can overcome the damping forces. 

Using equations (69)-(70) to eliminate the total pressure perturbation, the gas momentum 
equation (47) for short wavelength acoustic perturbations may be written 
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(71) 



The destabilizing force density is k(k-F)<5p times a numerical factor depending on opacity derivatives 
that may be positive or negative. For the purposes of discussion, let us assume for the moment that 
it is positive. Figure 4(a) illustrates the geometry of this force acting on a downward propagating 
acoustic wave. In this case k(k • F)bp is everywhere opposite to the velocity perturbation, and the 
wave is therefore damped. However, for an upward propagating wave, shown in Figure 4(b), the 
force is everywhere in phase with the velocity perturbation, and therefore amplifies the wave. If 
we had chosen the numerical factor depending on opacity derivatives to be negative, then it would 
have been the upward propagating wave that is damped and the downward propagating wave that 
is amplified. 
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It is also possible to consider the effects of the destabilizing force density in terms of whether 
it induces a net lag or lead between the total pressure and density perturbations in the sound wave. 
Once again, only one direction of propagation is destabilizing when viewed in this way. 

Examination of the growth rate in equation (57) reveals that it is just the combination of 
opacity derivatives P + (7 — l)0Tg that determines the action of the force, not the full factor 
[1 + P + (7 — l)0Tg] in equation (71). Pure Thomson scattering media are not unstable by 
this mechanism. To see why this is so, it is necessary to examine the equation of motion (71) 
more closely. Using the perturbed continuity equation (46) and regrouping, equation (71) may be 
rewritten as 
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(72) 



It is convenient to rewrite this equation in terms of the Lagrangian displacement £, where 

-iuj£. (73) 



5v= d~t 



Using the fact that to lowest order in k 1 , Sp = —ipV. ■ | from the continuity equation (46), we 
obtain 



dt 2 



= -kc*(k-£) 



— 1 



^k(k-F) (k.€)-^k€-Vp + g(k-€) 



+ O(k M\. 



(74) 



The destabilizing k-F term that would be present in Thomson scattering atmospheres competes 
with the equilibrium gravitational acceleration and gas pressure gradient in the last term in square 
brackets in equation (74). Figure 4 illustrates this by showing that the projection of g5p along the 
wave acts in the opposite direction to k • FkSp. Eliminating the gas pressure gradient using the 
hydrostatic equilibrium equation (21), the last term in square brackets in equation (74) may be 
written 

^k(k-F)(k-€)-k€-Vp + g(k-€) = ~k(kxF)-(kx^)-€x (kxg). (75) 

Because hydrodynamic acoustic waves are longitudinal to lowest order in k^ 1 , £ oc k, and the term 
involving the flux vanishes. In addition, the second term does no work as it is perpendicular to the 
displacement. Hence these terms contribute nothing to the growth rate - the destabilizing radiation 
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pressure force that would be present in a Thomson scattering medium is exactly cancelled due to 
the hydrostatic equilibrium of the unperturbed flow, at least for longitudinal waves. In the high-fc 
limit, the modes with density perturbations that can couple to the equilibrium radiation force are 
necessarily longitudinal. This is why gravity waves do not experience radiative driving in the short 
wavelength limit. 

Taking the scalar product of equation (74) with the wave vector k, we obtain 

d 2 „ \ ,99/, \ , 9/ , . up I 4E ( oj \ AEcuJj — l) 2 

^(k-C) = -fc 2 c 2 (k^) + ^ 2 (k.^)^|— + 



dt 2 g c I 3 \k 2 ' KFpto 

- l(k-F)[e p + ( 7 -i)e Tg ]|+o(fc)|€|. (76) 

Replacing d/dt with — iu>, we immediately recover equation (57) for the mode frequency and growth 
rate. 

We may repeat this procedure to better understand the instability in the limit where the gas 
and radiation exchange heat extremely rapidly. The gas and radiation temperatures are then locked 
together, so that 

5T„ = ST, = ST = —SE = T ( — — — ^ . (77) 

In this regime we must replace the separate gas and radiation energy equations with a single, 
combined energy equation which, when perturbed, gives 

-iu(-^- + SE] = --^-5E-ioj(-^+3c 2 }5p-pT5v-V(S R + S T ) 
V7-1 / 3k f /> \7-l J 

+ ik-F [(i + e p )^ + (e Tg + G Tr )^l • (78) 



Once again, equation (78) immediately implies that 8E ~ 0(k~ 1 )5p for short wavelength acoustic 
waves. Rapid radiative diffusion smooths out both gas and radiation temperature fluctuations in 
this case. Equation (77) then implies that 5p = c 2 5p to lowest order, i.e. that the acoustic waves are 
isothermal in this limit. Equation (78) can then be solved for the radiation pressure perturbation, 

f = " V {" + T ) (f) " + wJ k ■ "X 1 + W + °^ s »- < 79 > 

From equation (77), the gas pressure perturbation is given to the same order by 

Sp = cfSp - ^ (l + %) p(%) S P + <J* (H) (k ■ F)(l + Q P )Sp + 0(k- 2 )S P . (80) 



The first term on the right hand side of equation (79) and the second term on the right hand 
side of equation (80) are the damping terms due to radiative diffusion. The last terms in both 
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equations are the potentially destabilizing terms from the equilibrium radiation flux acting on 
density and opacity fluctuations. Note that the opacity derivative with respect to temperature, 
@Tg, has been lost because the gas temperature is locked to the uniform radiation temperature. 
There is also no buoyancy response in the pressure perturbations to this order, as expected because 
of the suppression of gravity waves in this limit. 

Using equations (79)- (80), the gas momentum equation (47) becomes 

d6v 



at 



+ 



k 2 c 



■*'-T k ( 1 + §)(T + ')(£)«' 

k (k ■ F) (l + |0 (1 + @ p ) bp + gbp + Oik-^Sp. 



(81) 



Once again, the destabilizing force density term is proportional to k(k • F)5p, and the geometry of 
Figure 4 still applies. Repeating the same steps as before, we obtain 
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which should be compared with equation (74). The destabilizing k • F term that would remain in 
a pure Thomson scattering atmosphere again competes with other forces in the last term in square 
brackets. Using the equilibrium LTE condition (20), the hydrostatic equilibrium equation (21), and 
the radiative diffusion equation (23), the equilibrium density gradient is given by 



3p \ kf . 



cfVlnp = g + ^l + i |j -P. 

The last term in square brackets in equation (82) is therefore 
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(84) 

Once again, these terms do not affect the growth rate for longitudinal waves. After taking the 
scalar product of equation (82) with k, we finally obtain 
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(85) 



from which we recover the growth rate (62). 
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4.1. Relation to the NAR Approximation 



Gautschy & Glatzel (1990) have shown that strange modes are well-described by a so-called 
non-adiabatic reversible (NAR) approximation. Mathematically, this approximation reduces to 
taking the scalar product of the perturbed radiative diffusion equation (32) with k, and then 
setting k • <5F = 0. One also ignores the radiation energy equation (29). In addition, if the gas and 
radiation temperatures are locked together by rapid emission and absorption of photons, the gas 
energy equation (28) is also ignored. The perturbed radiation pressure is then given by 
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(86) 



This equation recovers our high k expressions (70) and (79) for the perturbed radiation pressure, 
except that the NAR approximation cannot reproduce the radiation diffusion damping terms (as 
it neglects the radiation energy equation). Indeed, if we use equations (32), (70) and (79) to derive 
the perturbed radiative flux, we find 
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' (87) 
l + Q P )i + 0{k- l )8p if ST r = 8T g . 



When unstable driving (the second term on the right hand sides) strongly dominates radiative 
diffusion damping (the first term on the right hand sides) , the NAR approximation that k • 5F = 
is excellent. 

By not invoking the NAR approximation, however, we are able to deduce the threshold criteria 
listed in Table 1 required for the unstable driving to exceed the damping by radiative diffusion. 
At the same time, we can also calculate the actual luminosity perturbation associated with the 
unstable wave. 



5. MHD Instabilities in a Static Medium 

We now extend the analysis of the previous section by including magnetic stresses. It turns 
out that magnetic fields widen the domain of the acoustic wave instabilities to include even pure 
Thomson scattering media. This is because magnetic tension endows acoustic waves with a mixed 
longitudinal/transverse character, except in the special cases of propagation purely along or per- 
pendicular to the magnetic field. 

5.1. Short Wavelength Limit 

The dispersion relation (49) again factors easily in the short wavelength limit: 
= [u 2 - (k • v A ) 2 ] [u, 4 - ^k\v\ + c 2 ) + k 2 4(k ■ v A ) 2 ] 
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From left to right, the factors correspond to Alfven waves, fast and slow magnetosonic waves, the 
damped radiative diffusion mode, and a damped 5T g ^ 5T T mode. 5 The gravity waves of equation 
(54) have been lost because they are dominated by magnetic fields in this short wavelength limit, 
at least when k • va 7^ 0, i.e. for wave vectors that are not orthogonal to the equilibrium magnetic 
field. 

The first order corrections to the Alfven modes vanish, so that to = ±k • va + 0(k^ 1 ). Be- 
cause Alfven waves are purely transverse and thus do not involve density fluctuations in this short 
wavelength limit, they do not couple to the radiation physics. 

To the same order, the magnetosonic wave frequencies are given by 



= ±kv ph - 



2cv ph 



J ph 



(k • v A ) 



2v 2 



f AEv ph | 4£ccj a ( 7 -l) 2 ^ 



KFpVph 



(k-v A )(kx v A ) ■ (k x F) + ©(AT 1 ) 



' 2cv ph( 2 ^ h -v 2 A 



k-F) [e p + ( 7 -i)e rg ] 

(89) 



where > is the phase speed of the fast and slow magnetosonic waves, given by 



= i R + cl ± (v\ + elf - 4(k • v A fcl 



1/2 



(90) 



Equation (89) bears a strong resemblance to its hydrodynamic counterpart (57). The gas sound 
speed c g has been replaced by the phase velocity of the relevant magnetosonic wave v p h, and the 
"hydrodynamic" corrections to the mode frequency have been multiplied by a factor 



(k • va) 



zu ph 



> 0. 



(91) 



The last term in equation (89) is a new, destabilizing term that exists only in the presence of 
magnetic fields. As we discuss in more detail in section 6 below, its physical origin ultimately 
arises from the ability of magnetic stresses to support velocity perturbations that are not purely 
longitudinal, i.e. that are not parallel or antiparallel to the wave vector k. Note that this term 
vanishes for wave vectors that are either parallel or perpendicular to the equilibrium magnetic field. 
In each of these cases the compressive magnetosonic wave is longitudinal. 



5 In our previous paper (BS01), we found only seven nonzero frequency modes in our analysis of the perturbed 
radiation MHD equations. This is because we did not include heat exchange between the gas and radiation through 
absorption and emission, a process which turns out to be important in most applications (see section 7). Setting 
cj a = in equation (88) produces an additional zero frequency mode. 
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In the rapid heat exchange (u; a — > oo) limit, one of the eight modes is lost for the same reason 
as in the hydrodynamic case, and the remaining seven modes factor in the short wavelength limit 
as 
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Once again, instabilities only occur in the fast and slow magnetosonic waves, whose frequencies are 
given to first order by 
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(93) 



where v p h is given by equation (90) with c g replaced by the isothermal sound speed q. 



Equation (90) implies that (2-u 2 h — v \ — c 2 ) is positive for fast modes and negative for slow 
modes. Hence if the last "magnetic" destabilizing term in equation (89), or (93), dominates the 
hydrodynamic terms, the unstable fast and slow waves will have opposite vector phase velocities. 
It is also easy to see that in this case the slow mode will always grow faster than the fast mode. 
However, if hydrodynamic driving (the opacity derivative terms in eqs. [89] and [93]) dominates 
magnetic driving, then whichever of the two modes has the larger density perturbation for a given 
velocity perturbation will be the one that grows fastest. If v\ 3> Cg, the slow mode grows faster, 
whereas if v\ <C c 2 , the fast mode grows faster. 



5.2. Short Wavelength Limit With k • va = 0: Recovery of Gravity Waves 



If we consider wave vectors that are orthogonal to the equilibrium magnetic field, then the 
short wavelength limit recovers the fast magnetosonic modes, the radiative diffusion mode, and 
gives us a cubic equation describing coupled gravity and 5T g ^ 5T Y modes: 
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The reason we recover gravity waves in this limit is that magnetic tension does not exist when 
k • va = 0, and therefore no longer dominates buoyancy at short wavelengths. For weak thermal 
coupling between the gas and radiation, uo & — > 0, the roots of this dispersion relation correspond to 
magnetically altered gravity waves in the gas, 
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and a mode 
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In the rapid heat exchange limit, w a — > oo, the last mode is lost and the gravity waves depend only 
on the vertical density profile, 
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5.3. The Limit of Zero Gas Pressure: Photon Bubbles 



Just as in the hydrodynamic case, a number of authors (Arons 1992, Gammie 1998) have 
investigated the slow mode instability in the limit of zero gas pressure. Because the phase speed of 
the slow mode then vanishes, the instability then loses its obvious connection to this mode. 

Setting p = = c 2 in the dispersion relation (49), and then assuming uj oc as k — > oo, we 
find two modes with 
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(98) 



Just as in the hydrodynamic case [eq. (64)], equation (98) can be obtained from the thermally 
locked slow mode frequency in equation (93) by first squaring the frequency and then taking the 
limit of zero gas pressure. Again, we lose the damping terms in this limit, which determine the 
instability threshold for short wavelength modes. 

If we first assume negligible thermal coupling between the gas and radiation, and then take 
the limit of zero gas pressure, we obtain instead 
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(99) 



In a Thomson scattering medium, @ p = B^ g = 0, and equations (98)- (99) are identical to 
the photon bubble dispersion relation derived by Gammie [1998, his eq. (41)]. As noted both 
by Gammie (1998) and BS01, this growth rate only applies for wavelengths longer than the gas 
pressure scale height ~ c^/g. Unless this wavelength scale is optically thin, in which case our 
diffusion-based analysis breaks down, then the maximal growth rate occurs at wavelengths shorter 
than the gas pressure scale height, and is given by equations (89) or (93). 
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There has been some confusion as to how Gammie's (1998) work on photon bubbles relates 
to the original analysis performed by Arons (1992), who also assumed a pure Thomson scattering 
medium with zero gas pressure. Gammie's (1998) result can in fact be obtained as the short 
wavelength limit of Arons' (1992) dispersion relation [his eq. (33)], at least within the latter's 
assumed vertical magnetic field geometry. Arons' (1992) numerical growth rates all climb toward 
higher wavenumber, and must therefore eventually recover the k 1 ! 2 dependence of equation (98), 
provided the medium remains optically thick at these short wavelengths. Arons (1992) was mainly 
interested in long wavelengths where radiation diffusion was slow compared to the (radiation) 
acoustic wave period. Gammie (1998), on the other hand, was working in the rapid diffusion 
regime at shorter wavelengths. 

5.4. Numerical Results 

Table 2 summarizes the instability thresholds, growth rates, and characteristic wavenumbers 
for the magnetoacoustic wave instabilities in media with a variety of ratios of magnetic, gas ther- 
mal, and radiation energy densities. We have assumed a pure Thomson scattering equilibrium 
throughout this table, so that the hydrodynamic driving forces that we addressed in sections 3 and 
4 vanish, and only the magnetically-based driving remains. 

Figure 5 illustrates growth rates of unstable slow waves for small values of the thermal cou- 
pling frequency to t h for a radiation pressure dominated equilibrium, for which the asymptotic high 
wavenumber growth rate is accurately given by equation (89). Just as in the hydrodynamic case, 
the short wavelength growth rate declines with increasing thermal coupling. This continues until 
Uth exceeds g/c g , beyond which the gas and radiation temperatures are effectively locked and the 
waves reach unstable asymptotic growth rates given by equation (93), as shown in Figure 6. 

Just as in the hydrodynamic case, the unstable growth rate for high thermal coupling only 
extends over a finite range of wavenumbers. Above a critical cutoff wavenumber, the waves again 
become damped. After numerical experimentation, we deduced the physical dependencies of these 
cutoff wavenumbers, and they are listed in Table 2. Figure 7 illustrates the accuracy of these 
formulas in two specific cases, one with va > c g and one with v\ < c g . 

6. Physics of the MHD Wave Instabilities 

The magnetoacoustic wave instabilities are much more robust than their hydrodynamic coun- 
terparts, and in particular they survive unscathed even in pure Thomson scattering media. In order 
to see how this arises physically, we repeat the analysis of section 4 here to re-examine the com- 
peting forces acting on a parcel of gas in the wave. There are now two relevant waves that possess 
density fluctuations: the fast and slow magnetosonic waves. For reference, their polarizations arc 
shown in Figure 8. 
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Using equation (48) to eliminate the magnetic field perturbation, and equations (69)-(70) to 
eliminate the total pressure perturbation, the perturbed gas momentum equation (47) may be 
written in the high wavenumber limit as 
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Apart from the additional magnetic terms, this equation is of course identical to its hydrody- 
namic counterpart, equation (74), which we discussed above in section 4. There we found that the 
last term in square brackets did not affect the growth rate. This continues to be true for the piece 
of this term that involves the gravitational acceleration g, as this piece does no work. However the 
first piece that involves the radiative flux now can play a role because magnetosonic waves are, in 
general, not longitudinal due to the effects of magnetic tension. 

In the short wavelength limit, the fast and slow magnetosonic modes described by equation 
(100) are polarized in the plane of k and va, and their Lagrangian displacements may be written 

as 

€ = [v 2 ph k - (k • v A )v A ]V = e*/>, (101) 

where tp is some complex scalar amplitude. Taking the scalar product of equation (100) with the 
polarization vector e, we obtain 
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where, in the last term, we have used the fact that (k • e)(k x £) = (k • £)(k x e). The expression 
for the growth rate, equation (89), now follows immediately from equation (102) and the fact that, 
to lowest order in k~ l , 

~t = (103) 

As illustrated in Figure 8, the polarizations of the fast and slow waves for fixed k are orthogonal, 
and this normalized projection of £ onto k therefore always has opposite signs for the two modes. 
Hence the unstable fast and slow modes propagate in different directions if the magnetic driving 
term dominates over the hydrodynamic terms. 
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To understand the physics of the MHD driving further, consider again the relevant force density 
(p times the second to last term in equation [100]), which can be written as 
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(104) 



MHD driving results fundamentally from a breakdown in cancellation of the two force densities on 
the right hand side of this equation, whenever the waves are not longitudinal. The first, (Kp/c)k(k- 
F)<5p, is the driving force based on density, not opacity, fluctuations that we illustrated in Figure 
4. The second arises from a change in the radiation pressure along a fluid displacement, and is 
z(«F/Vc)k(£ ' F). If we ignore the radiation diffusion damping and opacity fluctuations in equation 
(70), these two force densities simply combine to give the negative gradient of the Lagrangian 
perturbation in radiation pressure AE/3, 
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For any purely longitudinal wave (£ and k parallel), the two force densities exactly cancel, and 
the Lagrangian radiation pressure perturbation then vanishes (again ignoring radiation diffusion 
damping and opacity fluctuations). For an MHD wave with fixed density perturbation bp and fixed 
k, magnetic tension can alter the direction of £ off of k, provided propagation is neither along or 
perpendicular to the field. Spatial gradients in the flow (which determine the radiative heat flow) 
are no longer parallel to fluid displacements. The second force density on the right hand side of 
equation (104) will now be greater or less than the first, depending on the type of MHD wave (slow 
or fast) and the propagation direction. Note that the two force densities also cancel for vertical 
propagation. 



The physics is very similar in the case where gas and radiation exchange heat rapidly. The 
perturbed gas momentum equation is then 
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Taking the scalar product of this equation with e, we obtain 
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With the help of equation (103) with c g replaced by the isothermal sound speed q, this now gives 
us the growth rate of equation (93). 

We have extended Hearn's (1972) hydrodynamic acoustic wave analysis of optically thin me- 
dia to the magnetosonic wave instabilities, and have found that the extra magnetic piece of the 
instability vanishes completely in this regime. Physically, this is easy to understand. In the limit of 
vanishing optical depth, the dynamics of the radiation field can be ignored, and the equilibrium flux 
F merely acts to globally reduce the equilibrium gravitational acceleration g to g c g- = g — KpF/c. 
Hence, apart from the hydrodynamic driving terms that are proportional to logarithmic derivatives 
of kf, the radiation flux enters the gas momentum equation only through a(x(kx g c g) term [the 
last term in equations (100) and (106)]. This term performs no work. Presumably, there will still 
be some level of instability in media with optical depths < 1, but an analysis of this situation is 
beyond the scope of the present paper. 



7. Astrophysical Applications 

We intend to investigate applications of the theory presented in this paper to particular astro- 
physical phenomena in future work. For now, we limit ourselves to a brief discussion of how these 
instabilities might manifest themselves in accretion disks and stars. 



7.1. Accretion Disks 

A number of authors have considered applications of strange modes (Glatzel & Mehren 1996, 
Mehren-Baehr & Glatzel 1999) and magnetic photon bubbles or magnetoacoustic waves (Gammie 
1998; BS01; Begelman 2001, 2002) to radiation pressure dominated accretion disks around black 
holes. 

The gas temperature is probably reasonably well-locked to the radiation temperature in stan- 
dard accretion disk models. For stellar mass black holes, the most important atomic absorption 
opacity is bremsstrahlung, for which the Planck mean opacity is given in cgs units by 

k p ~ 2.2 x W 24 gppT- 7/2 , (108) 

where 
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gp = -— I dvg s e- hv,k ^ (109) 

is an appropriate frequency-averaged Gaunt factor (Rybicki & Lightman 1979). Equation (108) 
gives a lower limit to the Planck mean opacity, which should also include bound-free and bound- 
bound contributions. These contributions are generally more important for the cooler disks thought 
to exist around supermassive black holes in active galactic nuclei. 
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By approximating gg to be nearly constant inside the frequency integrals, it is also straight- 
forward to show that for bremsstrahlung, 

1 d In kj 7r 2 . 

^isM;^ 1 ' 6 ' (110 > 

For concreteness, consider a standard Shakura &: Sunyaev (1973) disk model in which the 
anomalous stress is a times the total pressure. Using equations (30) and (66), the thermal coupling 
frequency due to bremsstrahlung is 
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Here L is the disk luminosity, -LecM * s the (Thomson scattering) Eddington luminosity, M is the 
black hole mass, rj = L/(Mc 2 ) is the disk radiative efficiency, r is the radius of the particular point 
of interest in the disk, and r g = GM/c 2 is the gravitational radius of the black hole. We have 
neglected all relativistic correction factors in the disk model. The particular numerical value of the 
thermal coupling frequency in equation (111) assumes vertically averaged, disk interior quantities. 

The corresponding thermal coupling frequency due to Compton scattering is 
^ s KTfx (Ml) K 2 x 10 . rad .-.„-. ( |L ) 1 ^ - . (1 12) 

Comparing equations (111) and (112), we see that Compton scattering generally dominates 
the thermal coupling of gas and radiation in the inner regions of the disk, while bremsstrahlung 
becomes more important further out, at least in the disk interior for this class of models. For 
hydrodynamic acoustic waves, the requirement that the gas and radiation temperatures be locked 
in the wave is that w t h exceed the sound wave frequency on scales of order the gas scale height 
c 2 /g, i.e. that 0JthC g /g ^ 1. For standard Shakura & Sunyaev (1973) models, 

3 / \ 9/2 
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and 
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Hence while bremsstrahlung may be unable to ensure tight thermal coupling between the gas and 
the radiation in the innermost radii, Compton scattering is easily adequate. We stress, however, 
that these estimates are very uncertain as they do not properly take into account the vertical 
stratification in the disk. In particular, we have assumed a gravitational acceleration appropriate at 
the disk scale height but a temperature appropriate to the disk midplane. Because Wth.CompCg/fl' oc 



T 3 / 2 , it could be that in regions near the disk photosphere, the gas temperature may not be able 
to track the radiation temperature. 



If the disk is treated hydro dynamically, then unstable acoustic wave growth rates are likely to 
be much slower than the dynamical frequency in the innermost parts of the disk because Thomson 
scattering dominates the Rosseland mean opacity. From Table 1, the maximum growth rate for 
hydrodynamic instabilities is roughly given by 



for a radiation pressure supported disk with K a ~ Kp. 

For example, Glatzel & Mehren (1996) found that a class of global strange mode instabilities 
was present at 25 km radii in accretion disk models around 1 M Q black holes with accretion rates of 
10~ 9 M Q yr _1 . Our estimate (115) then gives a growth rate ~ 5 x 10 — 3 S7, in remarkable agreement 
with their calculated value ~ 3 x 10 3 Q. Mehren-Baehr <fe Glatzel (1999) also note that this class 
of modes has growth rates that decline as the ratio of gas pressure to radiation pressure increases, 
which is also in agreement with our estimate (115). 6 In addition to these "hot" instabilities, Glatzel 
& Mehren (1996) and Mehren-Baehr & Glatzel (1999) found a class of unstable modes in colder 
regions of disk models where gas pressure was substantial. They found that these modes required 
the presence of a vertical opacity peak due to helium or hydrogen ionization, but that these modes 
also satisfied the NAR approximation and were therefore likely to be related to strange modes. 
We speculate that these modes are related to the local, WKB driving we find in the gas-pressure 
dominated case. Radiation pressure support is not essential for strange modes provided the gas 
and radiation temperatures are locked, as these authors implicitly assumed in their analysis. 

It is almost certainly true, however, that the presence of magnetic fields cannot be neglected 
in these flows, and these allow in principle much faster growth rates through the slow mode 
instability (Gammie 1998, BS01). From Table 2, the maximum growth rate in this case is ~ 
f2(c r /c g )min[l, (v^/cg)], which can be much faster than the orbital frequency in the radiation pres- 
sure dominated inner regions. The characteristic length scale of the instability is given by the 
reciprocal of the turnover wavenumber, which is the gas scale height. Again in the radiation pres- 
sure dominated inner regions, this is smaller than the disk scale height by the ratio of the gas to 
radiation pressure. It is far from clear how this instability might, if at all, manifest itself in the 
presence of MRI turbulence. The MRI acts on longer time scales (the orbital time), and possi- 
bly length scales as well, and this separation of scales between the two instabilities suggests that 
they might both be simultaneously present. Numerical simulations of vertically stratified, radiating 
shear flows will be required to address this question. 



6 Mehren-Baehr & Glatzel (1999) actually state that the growth rates are inversely proportional to [3 = 3p/E, but 
they do not present direct quantitative evidence for a dependence on cj?/Cg rather than c r /c g . 
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7.2. Stellar Envelopes on the Main Sequence 

Except possibly for short wavelength acoustic waves near the stellar photosphere, the gas and 
radiation temperatures are tightly coupled together in main sequence stellar envelopes. Thermal 
locking at the turnover wavenumber ~ g/cf requires u; t h ^ d/ci- For gas pressure supported 
envelopes, for example, this translates to the requirement that 

(±) (*L) 1 (HE*) ( P ^ ( 2 ^ ) _1 > l (ll 6 ) 

\L q )\M q ) V 4 / Vl0- 7 gcm-V ^minutes/ 1 ~ ' 1 j 

which, given realistic envelope opacities, is usually satisfied. 

From Table 1, if the gas and radiation temperatures are locked together, then the vertical 
radiative flux will destabilize short wavelength hydrodynamic acoustic waves in a radiation pressure 
dominated stellar envelope if FQ p > Ec{. It is convenient to express this in terms of the flux mean 
(Rosseland) optical depth rp of the particular envelope layer in question. Radiative diffusion implies 
that F ~ Ec/tf, so that short wavelength acoustic waves will be unstable in a radiation pressure 
dominated stellar envelope for all layers with optical depths satisfying 

TF < (-) @ P . (117) 



This can also be written as 



In a gas pressure supported stellar envelope, the instability criterion in Table 1 can be written 
in a similar fashion, 

"s(f)(fK < 119 > 

Now, p/E ~ L-^dd/L if gas pressure dominates radiation pressure, where L^dd is the Eddington 
luminosity appropriate for the Rosseland mean opacity kf- Hence the instability criterion (119) 
may be rewritten as 

This equation strongly suggests that even in gas pressure supported envelopes on the mid- to upper 
main sequence, acoustic waves may be unstable and strange modes might exist. 

We stress that we have only performed a local analysis of the stability of propagating waves in 
this paper. A demonstration that such waves are amplified by the background radiative flux does 
not necessarily imply that global compressive modes of the star will be unstable. This is because 
such modes can in general be viewed as superpositions of WKB waves propagating in different 
directions and, as we have seen, while one propagation direction may be driven, the other may very 
well be damped. Our instability criteria are therefore necessary, but not sufficient when it comes 
to global standing waves of the star. Note that the radiative driving becomes most important at 
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lower optical depths, and it is also important to point out that we have not included the effects of 
turbulent damping associated with outer convection zones, when they exist. 

Most main sequence stellar envelopes probably contain rather weak magnetic fields, and there- 
fore the hydrodynamic driving that we have discussed so far is probably of greatest importance. 
There are exceptions, however, where MHD effects might play a role. In some chemically peculiar 
A and B-stars, for example, highly ordered, roughly dipolar magnetic fields of order 1 kG have 
been measured (e.g. Landstreet 1992). The MHD driving that we discussed in sections 5 and 6 
above could produce instabilities in the upper layers of such stars, and magnetic effects will almost 
certainly alter the hydrodynamic driving. A particular case in point is that of roAp stars, which 
exhibit high order p-mode oscillations with low angular momentum quantum number £, largely 
confined to the poles of a dipole magnetic field (e.g. Kurtz 1990). Consider the possibility that 
hydrodynamic driving due to a large opacity variation dominates magnetic driving, but in the 
surface layers where the magnetic energy density dominates the gas pressure. From equation (90), 
the phase speed of slow modes in this case is ~ |k • va/^aIc;, assuming the gas and radiation 
temperatures are locked together. Equation (93) then implies a growth rate due to hydrodynamic 
driving which is proportional to |(k • VA)(k • F)OJ. In other words, slow waves will have maximal 
driving if they are vertical (i.e. radial) in the regions of the magnetic poles, where the field is also 
vertical. This might therefore be the basic driving mechanism behind roAp oscillations. We intend 
to explore this further in future work. 



8. Conclusions 

We have examined the conditions for local driving of a broad class of instabilities of propa- 
gating waves in optically thick media. These waves can be driven unstable by a sufficiently strong 
equilibrium radiative heat flux acting on density fluctuations in the wave. The central mathemati- 
cal results of this paper are the local dispersion relations: equation (49) for the full MHD case and 
equation (53) for hydrodynamics. Short wavelength expressions for the hydrodynamic damping 
and driving of acoustic waves are given in equations (57) and (62). Generalization of these for the 
full MHD case are given in equations (89) and (93). In addition to the acoustic waves, we have also 
briefly discussed short wavelength gravity waves in sections 3 and 5. 

There are essentially two types of local radiative driving of the acoustic wave instabilities. 
Hydrodynamic driving, which occurs even in the absence of magnetic fields, requires that the 
wave possess fluctuations in the flux mean (i.e. Rosseland mean) opacity. A medium with pure 
Thomson scattering opacity possesses no such local driving. The fastest growth rates occur for 
waves with vector phase velocities parallel (or possibly antiparallel) to the equilibrium flux F. 
Table 1 summarizes the properties of acoustic wave instabilities subject to hydrodynamic driving. 
It is this driving that is responsible for the existence of global strange mode oscillations in stars, 
which have normally been found in radiation pressure dominated stellar envelopes. We find that 
hydrodynamic driving can also produce fast growth rates in propagating acoustic waves in gas 
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pressure dominated envelopes, at least for sufficiently large radiative heat fluxes and provided that 
the thermal emission and absorption effectively lock the gas and radiation temperatures together. 

The presence of a large scale, equilibrium magnetic field expands the ability of an equilibrium 
radiative flux to drive acoustic waves unstable. In particular, MHD driving requires only that the 
wave possess density fluctuations, and even a medium with pure Thomson scattering opacity can 
have unstable magnetosonic waves. If MHD driving dominates over hydrodynamic driving, then 
the fastest growth rates occur for slow magnetosonic waves, and it is straightforward to show that 
these growth rates occur for phase velocities that are in the plane of the equilibrium flux F and 
magnetic field B. In general, the fastest growing waves have phase velocities that are at some 
angle to both F and B, and the MHD driving in fact vanishes for propagation along F, along 
B, or perpendicular to B. The characteristics of MHD driving in various limits are summarized 
in Table 2. MHD driving is responsible for "photon bubble" modes in accreting X-ray pulsars 
and accretion disks around black holes and neutron stars. It appears to require diffusive radiation 
transport: MHD driving vanishes in the optically thin limit, in contrast to hydrodynamic driving 
which survives unscathed (Hearn 1972). 

It it important to note that magnetic fields can also alter the behavior of the hydrodynamic 
driving that exists in the presence of opacity fluctuations, and this is summarized in equations (89) 
and (93). We suggest that this magnetically modified hydrodynamic driving may play a role in the 
excitation of observed p-mode oscillations in roAp stars. 

The physics of radiative driving may be more generic than just photon-matter interactions. In 
a coupled two-fluid system, provided one of the fluids dominates the inertia and the other provides 
rapid diffusive heat transport through which momentum exchange can occur, then unstable driving 
of density fluctuations of the sort we have discussed here may occur. One example might be diffusive 
neutrino transport in proto-neutron stars or hyper-Eddington accretion flows. Rapid streaming of 
cosmic rays in the interstellar medium can also drive acoustic waves unstable (Begelman <fc Zweibel 
1994), though it is unclear whether and how this might be related to the instabilities we have 
considered here. We hope to explore more detailed applications of radiatively driven instabilities 
to particular astrophysical phenomena in future work. 
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A. Appendix: Additional Mathematical Justification for the WKB Damping and 

Driving Terms 

Our WKB treatment presented in sections 3-6 is not rigorous, particularly as the interesting 
physics (the damping and driving effects on the underlying acoustic and magnetoacoustic waves) 
lies in a first order correction to the infinitely short wavelength limit. Because we invoked the 
WKB approximation on our separate perturbation equations before attempting to combine them 
into one, these corrections are vulnerable to alteration by other equilibrium vertical gradients. 

A more rigorous approach would be to combine all nine of our original linear perturbation 
equations together into one single ordinary differential equation in z, and then apply WKB tech- 
niques on that equation alone. Unfortunately, we have not discovered a way of doing this in general. 
However, there are two particular cases where this can be done, and we present them here. The 
results of these cases fully agree with the asymptotic WKB results that we obtained in sections 3-6. 



A.l. Case 1: Vertically Propagating Hydrodynamic Acoustic Waves With Equal 

Gas and Radiation Temperatures 

Assume that the gas and radiation exchange energy quickly enough to guarantee that T r = 
T g = T, and neglect the effects of magnetic fields. Consider perturbations that have no variation 
in the horizontal direction, i.e. waves which propagate (up or down) in the vertical (z) direction. 
The linear perturbation equations for these waves are 

d5p . dp dSv n , . 

d5v d6p r IdSE lk , 

p -at = -te- 9Sp -s-»r> (A2) 

rr^ r ( du dE \ ( 4 85F 

188E k-fp KpF KpF 5T 

together with 5u = 5p/{^ - 1), 5p/p = 5p/p + 5T/T, and SE/E = A5T/T. Note that in this case 
of vertical propagation, 5v and <5F only have nonzero components in the vertical direction. 

These equations are still too complicated to combine into a single wave equation, and we must 
approximate the thermodynamics somewhat. The perturbed energy equation (A3) can be written 
just in terms of 5T, 5F, and 5v as 

If p ir \ fdST . dT\ ( \\dbv d5F 



T \7 — 1 / V dt dz ) V 3 ) dz dz 

Now, at very short wavelengths (k — > oo), we expect rapid radiative diffusion to smooth out 
temperature fluctuations, so that ST = 0(k~ 1 )Sv is very small. Hence the dominant terms in this 
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equation at short wavelengths are the last term on the left hand side and the term on the right 
hand side. 7 At the same level of approximation, we may therefore rewrite this equation as 



d_ 

dz 



p+-E)Sv + SF 



0, 



which can be immediately integrated to give 



p + -E ) Sv + SF = 0. 



(A6) 



(A7) 



Similarly, we also neglect all ST terms in the momentum and diffusion equations (A2) and (A4) 
that are not gradients to give 

dSv v dT „ pdSp pdST „ AE dST 



dt 



p dT pdbp pd5T 
"pT~dz~ 9 ~ ~p~dz~ ~T~dz~~ 9P ~%f~dz~' 



and 



respectively. 



4E d5T 
^T~dz~ 



w S F-^(i + e p) 5 P , 



(A8) 
(A9) 



With these approximations, equations (Al), (A7), (A8), and (A9) can now be combined into 
a single equation. Without further loss of generality, we assume all perturbation variables have 
a time-dependence of the form exp(— iujt). Then after some algebra, we obtain a single ordinary 
differential equation in Sv: 

d 2 5v 







dz 2 



+ 



dp 



dz 



?yp 



AE 



Kp P F 



d5v 
dz 



+ 



o lUJKYp 

OJ p-\ 



3p kfF dp K-pF 



@o 1 + 



3p\ dp p f dp 



AE c dz c " \ AE J dz p 2 \dz 
It is convenient to define a new perturbation variable Sij) by 

Sip = p x l 2 bv. 

Then equation (A10) becomes 



+ 



p d 2 p 
p dz 2 



AE 



Sv. 



1 + ^ 
AE 



(A10) 



(All) 



P- 



cPStp 
' ~dz^ 



- 1 + 



3p 
AE 



e, 



KppF dSip 
c dz 



+ 



AE\ 



3 J 



3p 



AE 



1 dp K-FpF 



2p dz 



1 + ^ 

AE 



+ _ dp 
Ap \dz 



ld 2 p 3p kfF dp KpF t 
Tdz 2 ~ ~AE~ir~dz 



AE dz 



p I dp\ p d? p 



dz 



+ 



•dz 2 



8il>. 



(A12) 



7 This approximation was also used by Begelman (2001) in his treatment of nonlinear, radiatively driven MHD 
waves, cf. his equation (7). 
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If we now employ the WKB ansatz that 5ip oc exp(ifez), this equation gives the dispersion relation 

+ O(k ), (A13) 



or 



LU 



LU 



,2 2, - K F f, . 3p 



4E\ 



1 + 



3p 
4E 



V + 



AE 



+ o(k- 1 ), 



(A14) 



in exact agreement with equation (62) for vertically propagating waves. 



An astute reader familiar with WKB techniques may notice that we cheated slightly The 
formal WKB treatment first involves eliminating the first order derivative term, thereby forcing the 
differential equation to resemble the harmonic oscillator equation. We can do this by transforming 
from 8v to Sift, defined by 



Sip = p l l 2 5v exp 



e, 



KppF 

lap 



(A15) 



When we employ the WKB ansatz to the differential equation that then results, 5ip oc exp(ifcz), we 
obtain 



±kc, — i — 
2c 



1 + lE 



P 



(A16) 



This appears to be an acoustic wave which is merely damped by radiative diffusion, and we have 
apparently lost our radiative driving term. However, one must be very careful in the physical 
interpretation of the wave vector k. It is related to the wave vector k we employed above by 



k = k + ie, 



KppF 

lap 



AE 



(A17) 



Substituting this expression into equation (A16), we recover the dispersion relation (A14), including 
the radiative driving term. In order to interpret this ambiguous behavior, note that the velocity 
perturbation is 



Sv oc p 1 I 2 ex.p(—iujt + ikz) oc p exp 



— iut + ikz + / G 



KFpF 

2cp 



(A18) 



For concreteness, assume Q p > and consider upward propagating waves. We may choose k to be 
real and large (the interpretation we have adopted throughout this paper), in which case equation 
(A14) is the correct dispersion relation. If radiative driving dominates damping, then the wave will 
grow in time. Alternatively, we may choose k to be real and large, in which case equation (A16) is 
the relevant dispersion relation. The wave will tend to damp with time, but if the same quantitative 
requirement that the driving term dominates the damping term in equation (A14) still holds, then 
equation (A18) implies that it will grow exponentially as it propagates spatially. Either way, we 
have instability: the wave is driven to larger and larger amplitude by the background radiative 
flux. 
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Finally, we note that the original partial differential equations (Al), (A7), (A8), and (A9) can 
be combined into the following energy conservation equation for the waves, 



1+^ 

4E 



K F F 



Qp5p5v 



4E\ KFP _ 6v2 



(A19) 



This equation leaves no ambiguity in interpretation. The two terms on the right hand side represent 
radiative driving and damping by diffusion, respectively, and are consistent with the dispersion 
relation (A14). Note that for Q p > 0, radiative driving requires SpSv > 0, i.e. the density and 
velocity perturbations must be in phase, which occurs only for upward propagating waves. This is 
also fully consistent with the physics we described in section 4 and illustrated in figure 4. 



A. 2. Case 2: Infinitely Strong Horizontal Field, Equal Gas and Radiation 

Temperatures 

Assume again that the gas and radiation exchange energy quickly so that they have the same 
temperature, and consider a case where there is a horizontal equilibrium magnetic field that is 
so strong (v^/cg — > oo) that fluid elements are constrained to move horizontally. In a Cartesian 
coordinate system (x, y, z), with x along the field and z in the upward vertical direction, the linear 
perturbation equations in this case may be written as 

d5p d5v n , . . 

ir + "*= ' (A20 > 

dSv d5p kfp f k \ 

pi*=-i£ + ^ SF - < A21 > 

_({„ + j£) + ^„ + _ E j _ = _V ■ ,5F, (A22) 

if = -^.. im 

185E K F /9 rr , kfF „ „ . . KppF 5T 

3^7 = —r SF > ^ )5p ~ ~ir eT T' (A25) 

together with Su = 5p/(-y - 1), Sp/p = 5p/p + 5T/T, and 8E/E = A5T/T. The magnetic field has 
vanished entirely from these equations, as it is considered to be so strong that fluid motions do not 
bend the field lines. As a result, Alfven and fast magnetosonic modes have been eliminated from 
this system of equations. Note that, in contrast to the previous case we considered, we have not 
made any approximations concerning the rate of radiative diffusion. 

Without loss of generality, we assume that all perturbation variables have an (x,y,t) depen- 
dence proportional to exp[i(k x x + k y y — ut)\. Then with some algebraic work, all perturbation 
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variables can be eliminated in favor of ST, giving us a single, second order ordinary differential 
equation in z, viz. 







4Ec d 2 



+ 



3kfp dz 2 \ T 
P 



' k 2 x F(l + Q p )(p + 4E/3) 
poj 2 - pk 2 

ik 2 x oj{p + AE/3) 2 



10J 



7-1 



+ F(Q T - 4) + — 

dz 



d ( AEc\ 



\3k f p) 



d_ f5T 

dz 



AEc 



(poj 2 - pk 2 ) 3n F p 



(k 2 + k. 



5T 
~T' 



(A26) 



(We have assumed for simplicity that Q p and Qt are constants.) Now, we expect that slow mag- 
netosonic modes will have a dispersion relation of the form oj 2 = k 2 c 2 = k 2 p/p plus higher order 
corrections. Hence adopt the WKB ansatz that 5T/T oc exp(ik z z) with large k x , k y , and k z ; lo oc k; 
and to 2 — k 2 p/p oc k. Then equation (A26) gives the dispersion relation 



Z 7 _ _ ' 

oj = k x c { + % 



2 2, -kx K FP 



k 2 



or 



oj = ±|A; x |ci 



2cc\ 



1 + 



3p 

1e 



AE 



k z F(i + e p ) 



4E\ oj 

" + T)- Pi 



P + f)«^(H9,, 



0{k c 



+ o(k- r : 



(A27) 



(A28) 



which agrees completely with equation (93) in the va — ► oo limit. This confirms that our crude 
WKB analysis presented in the main body of the paper also appears to work in the MHD case, 
and is robust with regard to the radiation physics. 



B. Appendix: Radiative Diffusion and Gas/Radiation Heat Exchange Effects on the 

Magnetorotational Instability 

In our first paper (BS01), in addition to pursuing an initial exploration of radiatively driven 
instabilities in magnetoacoustic waves, we also examined how radiative diffusion modifies the mag- 
netorotational instability (MRI; Balbus & Hawley 1991). Our analysis in that paper assumed pure 
Thomson scattering opacity, with no heat exchange between the gas and radiation. Here we wish 
to briefly consider the behavior of the MRI under the more general thermodynamic assumptions 
used throughout this paper. 

We restrict consideration to axisymmetric perturbations on axisymmetric equilibria rotating 
with angular velocity £l(R), where R is the radial distance from the rotation axis. Apart from the 
radial dependence of f2, which is the source of free energy for the MRI, we completely neglect all 
other equilibrium gradients in the flow. In order to avoid time-dependent azimuthal fields in the 
equilibrium that result from the shear, we assume that the radial component of the magnetic field 
is zero in the equilibrium. 
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Using cylindrical polar coordinates (R, <fi, z), the analysis proceeds very much as in BS01, with 
a resulting dispersion relation 



D ms D B u + k 2 z v 2 Az v A(j) I k z R 



dfl 2 
~dR 



i 2-2 

k to 



Here 



~2 _ 2 t 2 2 

U =u -k z v Az , 
D nis ^u> 2 -k 2 z (C 2 + v 2 A ,). 
k 2 

„ _ ~4 2-2 j n 2i 2 2 
£>BH = TfW - K U -Ail k z V Az , 
K z 

and C 2 = Ac 2 + Be 2 + C in the limit of zero stratification. From equations (41)-(44), this is 



(Bl) 

(B2) 
(B3) 

(B4) 



«-{(" + &)[" + 

icfc 2 \ o 
w ( w + , c 2 + 



4i(7 - l)uj a E 



P 



9( 7 - 



/ cfe 2 \ 



1 + 



4( 7 - 1)E 
P 



(B5) 



Here c s i is the total sound speed in the fluid, defined by 



2 _ 16(7 - l)^ 2 + 60(7 - l)Ep + 9jp 2 _T l (p + £7/3) 

c sl = 



9[p + 4( 7 - l)£]p 



(B6) 



where Ti is the first generalized adiabatic exponent commonly used for matter and radiation (Chan- 
drasekhar 1967). 

Equation (Bl) should be compared to equation (39) of BS01, which it closely resembles. Not 
surprisingly, altering the thermodynamics has produced only one difference: a change in the effective 
sound speed to C s . The basic conclusions of BS01 therefore continue to hold, with only minor 
quantitative changes: 

If the equilibrium azimuthal field is zero (v A< f, = 0) , then the &r = MRI mode (the "channel 
solution" ) is completely unaltered by radiation effects. 

If radiative diffusion is unimportant on the scale of the MRI critical wavelength, [c£1/(3kfPv Az ) <C 
1], then the MRI will be unaltered provided the azimuthal field energy density is smaller than the 
total thermal energy density in the fluid. To be more precise, the MRI will be unaltered pro- 
vided v A( p <C c 2 + c 2 if the gas and radiation are not thermally locked together on an orbital time, 
u a [l + 4(7 — l)E/p] < Q, or v A( p <C c 2 x otherwise. In either case, the essential nature of the MRI 
is unchanged provided the magnetic field is subthermal. 
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If radiative diffusion is rapid on the scale of an MRI critical wavelength, then radiation pressure 
no longer helps enforce incompressibility in the MRI. In a radiation pressure dominated medium 
where the magnetic field is largely in the azimuthal direction, even subthermal (with respect to the 
radiation) fields can then reduce the growth rate of axisymmetric MRI modes. Unless v 2 ^ <C c 2 
if gas and radiation are not thermally locked, or <C cf otherwise, then the growth rate of the 
MRI will be reduced. 
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Table 1. Order of Magnitude Conditions for Radiation Hydrodynamic Acoustic Wave 

Instabilities? 



Thermal 


Pressure 


Instability 


Asymptotic 


Turnover 


Cutoff 


Regime 


Support 


Criterion 


Growth Rate 


Wavcnumbcr b 


Wavenumbcr 



u; fc >u, th (<5T r ^5T g ) £»p F0 > E x max [ Cg , (^) (£)] (i) 

Wfc>WUl («r r ^«r g ) f«p F0>Bxmax[ Cg ,(^)(Q] ©(f)(f) e(^)(f) 

Wfc < u, th (5T r = <5T g ) F»p F© p >£ Ci e„(*) "(f) (^|! + 0pl 

u> k < wu, (5T r = «5T g ) F«p Fe P > pCi e„(*) "(f) (^|! + 0pl 



a Here is simply [fc| times the phase speed of the wave (c g or c- u depending on the thermal regime, for the hydrodynamic 
acoustic waves in this table), u> t ^ is defined by equation (66), u a is defined by equation (30), © = |© p + (7 — 1)0tv|, and © p and 
©T g are defined in equation (45). All the results in this table are valid only in the regime where the photon diffusion time across 
a wavelength is faster than the wave period, but yet the waves are still optically thick so that photon diffusion still applies. This 
requires ftpp(c g /c) < k < Kpp. 

b The minimum wavenumber at which the asymptotic growth rate becomes valid. (The growth rate typically tends to zero as fc 1 / 2 
as k — > below this wavenumbcr.) In some cases, the turnover wavenumber may be so small that the WKB approximation is violated. 
This is an issue primarily for gas pressure dominated equilibria (p > E), in which case our analysis is valid only for fc 2> g/c|. In 
addition, if © or |0 P | are less than unity, as can happen in Thomson scattering dominated media, the turnover wavenumber is less 
than the reciprocal of the gas scale height, and the waves can no longer be purely acoustic in nature at this low wavenumber. 



Table 2. Order of Magnitude Conditions for the Radiation Magnetoacoustic Wave Instabilities. 3 



Mode 


Thermal 


Magnetic 


Pressure 


Instability 


Asymptotic 


Turnover 


Cutoff 




Regime 


Pressure 


Support 
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Growth Rate (7^) 


Wavenumber 


Wavenumber 
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SLOW 
SLOW 
SLOW 
FAST 
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FAST 
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SLOW 
SLOW 
SLOW 
SLOW 
FAST 
FAST 
FAST 
FAST 



> w t h 

> w t h 

^k > W t h 

^k > W t h 

^k > w th 

^fc > w th 

^fc > w th 

^fc < W t h 

^fc < Wth 

^fc < W t h 

^fc < W th 

<^fc < W t h 

^fc < w th 

^fc < W t h 

<+>k < Wth 



S 2 /8tt >> p 
B 2 /8vr << p 
B 2 /8tt >> p 
B 2 /8vr << p 
■B 2 /8vr » p 
B 2 /8vr << p 
B 2 /8vr » p 
B 2 /8tt « p 
■B 2 /8vr » p 
B 2 /8tt << p 
B 2 /8tt » p 
B 2 /8vr << p 
B 2 /8vr >> p 
■B 2 /8vr « p 
S 2 /8tt » p 
S 2 /8tt << p 



£»p F>£xmax[ Cg ,(^) (£)] 

£»P F>(|) S xmax[, A ,(^)(4)] 

£«p F>Exmax[ Cg ,(^) (Q] 

£«P F>(|) £ Xmax[, A) (^)(4)] 

£»P F>Exmax[ % ,(^)(|i)] 

£»P F>(|) £ xmax[c g ,(^)g)] 

E«p f>£xmax[v A ,(^) (£)] 

£«P F>(|) Sxm ax[ Cg ,(^)g)] 

E»p F>Eci 

E»p F >^<\ev a 

E«p F>pc { 

E«P F>(^pv A 

E»p F>Ev A 

E»p F >(i\ Eci 

E«p F>pv A 

E«P F>(^-)p Ci 



{9va)/c% 

(f)*A* 

(f)(^A)/0 2 
5/*>A 

Cfftft)/cg 

(f)»/"A 

(f)(^i)/4 
(g-"A)/c? 

ff/ci 
(3^a)/c 2 

Cfft£)/c? 

S/*>A 
Cfft£)/c? 
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9/4 



(f)*M 
(f)*M 

f/^A 

9/c 2 



(f) 



A 
9/4 



OO 
CO 
OO 

00 

OO 
CO 



9/cf 


(a; th7d ) 1/2 /ci 


9/cf 


c? Kh7d) 1/2 hi 


9/cf 


Kh7d) 1/2 /ci 


9/cf 


c 2 Kh7d) 1/2 M 


9/vi 


(^th7d) 1/2 /f A 




Kh7d) 1/2 /ci 


9h\ 


(^th7d) 1/2 /fA 


|W 2 


(wth7d) 1/2 /ci 



a Once again, is |fc| times the phase speed of the wave. For simplicity, we have assumed that the flux mean opacity ftp is independent of density and 
temperature, so that the hydrodynamic driving terms of the instability vanish. This could be the case, for example, if Thomson scattering is the dominant 
form of momentum transfer between the gas and radiation. 
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Fig. 1. — Scaled growth rates of hydrodynamic acoustic wave instabilities as a function of scaled 
wavenumber, for different values of the thermal coupling frequency uJtt- The magnetic field was 
set to zero, and other equilibrium parameters were F = 300/9Cg and c r /c g = 10. Because of the 
dominant radiation pressure, we neglected gas pressure gradients and chose a flux mean opacity 
that satisfied KpF/c = g. We also chose the flux mean opacity to have derivatives appropriate for a 
Kramers type law: © p = 1 and ©T g = —3.5. The first plot (a) shows growth rates for the downward 
propagating unstable two-temperature waves that exist at high wavenumbers. In addition, there 
are upward propagating one-temperature waves that are also unstable and exist only up to a finite 
cutoff wavenumber. The second plot (b) is a blow up of the left hand region of plot (a) to show 
these modes more clearly. In both plots, the waves are assumed to be propagating purely vertically, 
i.e. the wave vector k is either parallel or antiparallel to the flux F. Such mode directions have 
maximal growth rates. 
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Fig. 2. — Growth rates (a) and phase speeds (b) for the hydrodynamic instability at different, 
high levels of the scaled thermal coupling parameter CgUth/g, which labels each curve. Except for 
uj t h, equilibrium parameters are the same as in Figure 1. The horizontal dashed line in figure (a) 
indicates the asymptotic growth rate predicted by equation (62), while the vertical dashed lines 
indicate the cutoff wavenumber (o>th<7|l + @ p \/cf ) 1//2 . The upper and lower horizontal dashed lines 
in figure (b) indicate the adiabatic and isothermal [q = (3/5) 1 / 2 c g ~ 0.77c g ] gas sounds speeds, 
respectively. 
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Fig. 3. — Growth rates for the one-temperature hydrodynaniic instability. Equilibrium parameters 
were chosen to be appropriate for an intermediate mass main sequence stellar envelope. The equi- 
librium gas pressure is much larger than the equilibrium radiation pressure so that L/L^d — 10~ 2 
and F z = WpCy. The flux mean opacity was chosen to be kf = lcm 2 g _1 with 0^ = 1 and 
@Tg = —3.5. The wave vector was chosen to be purely vertical. Notice that the growth rates are 
still dynamical even though the equilibrium radiation pressure is much smaller than the gas pres- 
sure. The tight thermal coupling between the gas and radiation elevates the dynamical significance 
of the radiation pressure perturbations, allowing for relatively large growth rates. The downward 
propagating instability at low k works off of opacity perturbations arising from temperature fluctu- 
ations, while the upward propagating instability at high k is consistent with the one-temperature 
results of Section 3.2. 
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Fig. 4. — Geometry of radiation pressure force due to density and opacity fluctuations for (a) 
downward and (b) upward propagating hydrodynamic acoustic waves. Depending on whether the 
multiplier of the radiation pressure force is positive or negative, one direction is damped and the 
other is unstable, provided the radiation pressure force is large enough to overcome other forms of 
damping. 
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Fig. 5. — Growth rates for the two-temperature slow MHD instability. Equilibrium parameters 
are the same as Figure 1 except for Q p = ®Tg = 0, k z /k = 0.556, and a vertical equilibrium 
magnetic field with va = 5c g . Again, increasing the thermal coupling between the gas and photons 
leads to damping in this two temperature regime. The fast mode instability does not appear since 
driving from the equilibrium flux does not in this case overcome the damping produced by radiative 
diffusion. 




Fig. 6. — One temperature slow MHD instability for different, high levels of thermal coupling. 
Apart from uj t h, parameter values are the same as Figure 5. 



- 53 - 




-54- 




12 3 4 

log[kcVg] 



Fig. 7. — One-temperature slow MHD instability with (a) the same parameters as for Figure 6 
except for different values of Wth> and (b) the same parameters except that a weak equilibrium mag- 
netic field was chosen with va = c g /5. The horizontal dashed lines indicate the asymptotic growth 
rate predicted by equation (93). The vertical dashed lines indicate the predictions of our analytic 
expressions for the cutoff wavenumber in Table 2: (a) (uth9 '/ 'q) 1 / 2 / ' c\ and (b) Ciiuj^gjv^) 1 ! 2 jv\. 



-55- 





Fig. 8. — Polarizations of the fast and slow magnetosonic modes. All vectors shown lie in the plane 
of the page. 



